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Abstract. - As a model of decohering environment, we show that quantum chaotic system be- 
have equivalently as many-body system. An approximate formula for time-evolution of the reduced 
density matrix of a system interacting with a quantum chaotic environment is derived. This the- 
oretical formulation is substantiated by numerical study of decoherence of two qubits interacting 
with a quantum chaotic environment modeled by chaotic kicked top. Like many-body model of 
environment, quantum chaotic system is efficient decoherer, and it can generate entanglement 
between the two qubits which have no direct interaction. 
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Interaction of a quantum system with environment cre- 
ates correlations between the states of the system and of 
the environment. These correlations destroy the super- 
position of the system states - a phenomenon known as 
decoherence [1]. This phenomenon is believed to be re- 
sponsible for quantum to classical transition. Decoherence 
is also a major obstacle for designing quantum computa- 
tional and informational protocol [2]. Therefore, a deeper 
understanding of the phenomenon is required to address 
the fundamental question like quantum-classical transition 
and, to develop quantum computational protocol. 

In general, the environment is modeled by many-body 
system, e.g. infinitely many harmonic oscillators in ther- 
mal equilibrium (Feynman- Vernon or Caldeira-Leggett 
model) [3], spin-boson model [4], chaotic spin-chain [5], 
etc. In another approach, random matrix model of the 
environment is used [6,7]. Random matrix theory has well 
known connections with quantum chaotic systems. Hence, 
some studies have concentrated on the possibility of hav- 
ing quantum dissipation and decoherence due to the inter- 
action with chaotic degrees of freedom [8] . Recently, as a 
model of decohering environment, single particle quantum 
chaotic system has been considered [9] . This paper shows 
that the kicked rotator, a well studied model of chaotic 
system, can reproduce the decohering effects of a many- 
body environment. In comparison to complex many-body 
system, this simple deterministic system is very convenient 
for numerical as well as analytical studies of decoherence. 
Hence, single particle quantum chaotic system warrants a 
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special attention as a model of decohering environment. 

In this Letter, we establish direct equivalence of sin- 
gle particle quantum chaotic environment and Caldeira- 
Leggett type model of many-body environment by pro- 
viding a rigorous but straightforward treatment. We keep 
our results vis-a-vis a recent study in which decoherence 
in a quantum system is investigated under the influence of 
a collection of harmonic oscillators environment [10]. Our 
derivation first assumes weak interaction between the sys- 
tem and the environment. Using the interaction strength 
as a small parameter, we perform perturbative theory cal- 
culation. By exponentiating the perturbative expansion, 
we get an approximate formula for the non-perturbative 
strong interaction effect of the environment on the system. 
The approximate formula is then justified by numerical 
evidences. In numerics, we study decoherence of two non- 
interacting qubits which are individually interacting with 
a common quantum chaotic environment. We use chaotic 
kicked top, a very well studied model of quantum chaotic 
system [11], as the environment. 

Most general form of the Hamiltonian of a system S, in- 
teracting with an environment E, is H = Hg + He + H$e, 
where H$ and He are the Hamiltonians of the system 
and the environment, respectively, and H$e is the system- 
environment coupling Hamiltonian. We assume through- 
out this Letter that the decoherence time is much smaller 
than the system characteristic time. Hence we can ne- 
glect any dynamics of the isolated system and can discard 
Hg. We consider kicked quantum chaotic system as a 
model of the environment, so the general form of the time- 
dependent system-environment Hamiltonian is : H(t) = 
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Is ® H E (t) + Hj(t) where £T B (t) = H x + H 2 6{t - n) 
and Hi{t) = aVs <S> Ve; Vs and Ve are the coupling 
agents of the system and the environment, respectively. 
Parameter a determines the strength of the interaction. 
We now assume that V E commutes with H\. This is 
a very natural assumption as it separates environment 
dynamics from the system-environment interaction. The 
corresponding time-evolution operator U in between two 
consecutive kicks is: U = UiUo = Ui(Is ® Ue), where 
the coupling part Ui — exp(—iaVs <8> Ve), Is is a unit 
matrix which indicates the absence of any dynamics in 
the system, and U E — exp(— iH\) exp(— ii^)- Further, 
we assume that the initial joint state of the system- 
environment is an unentangled pure state of the form 
|^se(0)) = IV's(O)) ® IV'e(O)). We measure decoherence 
of the system by its loss of purity P{n) = Trs[ps(n) 2 } 
which varies from 1, for pure state, to l/iVg for completely 
mixed state, N$ being the Hilbert space dimension of S. 
And ps(n) = ^e\^ SE(n)} SE(n)\ is the system reduced 
density matrix (RDM), where \^se(ti)) is the joint state 
of the system and the environment at time n. 

We are interested to study entanglement between the 
two non-interacting qubits due to their interaction with 
a common chaotic environment. Following Wootters, the 
entanglement is measured by computing the concurrence 
C[ps{n)] — max(A, 0), where A = \f\[ — ^f\^ — VA3 — 
and Aj's are the eigenvalues of the matrix R(n) = 

Ps(n)ps(n); ps{n) = {<J y ® a y )p s {n)* {a y <g> a y ), ps(n)* 
being the complex conjugation of ps{n) in the computa- 
tional basis {|00), |01), |10), 1 1 1) } [12]. The concurrence 
varies from 0, for separable state, to 1 for maximally en- 
tangled state. 

Our perturbation theory approach is reminiscent of 
the method followed by Tanaka et. al. in the context 
of entanglement production between two coupled chaotic 
systems [13]. We define the interaction picture of the 
joint density matrix ]Ts~E~(n) = Uo(ny psE(n)Uo(n), and 
of any arbitrary operator 0(n) = Uo(n)'OUo(n), where 
Uo(n) = (Is®U e)™, and 0(n) represents the free evolution 
of O. Thus, the time evolution of psE(n) is determined 
by the mapping : psE(n) = Ui(n)psE {n — l)Ui(ny . The 
perturbative expansion for Ui(n) is : 

Ui{n) = 1 - iaV(n) - ^a 2 V{n) 2 + 0(a 3 ), (1) 

where V(n) = Vs(n) <8> Ve(ti), Vs(n) and V E (n) are as 
usual free evolutions of Vs an d V E , respectively. We 
introduce the eigenvalues and the eigenvectors of Vs as 
Vs\s) — \ a \s). By tracing out the environment, we obtain 
perturbative expansion of RDM of the system S at time 
n in the eigenbasis of Vs upto 0(a 2 ) as : 
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Here, (•) = Ty{pe(0) •} is an average over the environ- 
ment. We extend the range of validity of the perturbative 
expansion by exponentiating Eq. (2) 

[ps(n)] ss , ^ exp[- ia(X s - \ s >)g(n) - a 2 (\ s - A s /) 2 /(«) 
-ta 2 (\ 2 s -X 2 s ,)cf>(n)][ps(0)} ss ,, (4) 

and our numerics suggest that this expression works well 
even in strong coupling regime. The same approxima- 
tion was also taken in random matrix study of quantum 
fidelity decay [14], which later justified by supersymme- 
try calculation [15]. This approximation is also verified in 
experimental study of fidelity decay in quantum chaotic 
system [16]. 

The most important term in Eq. (4) is the function 
f(n), which is responsible for the decay of the off-diagonal 
terms of the RDM, leading to the loss of coherence of the 
system. Hence, we call the function f(n) as "decoher- 
ence junction". We get rid of the term g(n) by redefining 
Ve(/) = V E (l) - (Ve(/)}. Then Eq. (4) will be exactly 
identical to the expression obtained for the environment 
consisting of a collection of harmonic oscillators [17] . How- 
ever, it does not prove the equivalence of both the expres- 
sions. In order to do so, we have to show that the de- 
coherence function /(n) is equivalent to the decoherence 
function of Ref. [10]. 

In the following, we closely look at the properties of the 
decoherence function f(n). From the definition of f(n) 
given in Eq. (3), and after doing some algebra, we get: 
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where Ce(1, I') is the correlation function of the uncoupled 
chaotic environment given by 

Ce(U') - (Ve(Z)Ve(O) - (Ve(0>(Ve(Z')>- (6) 

Here we assume that the environment is strongly chaotic, 
and the phase space of its underlying classical system is 
bounded. Following Ref. [13], we assume further phe- 
nomenological properties of C E (l 7 l') ■ (1) In strongly 
chaotic regime, the distribution function quickly becomes 
uniform in the phase space, hence Ce(1,1) — (Ve(/) 2 ) — 
(Ve(0) 2 becomes almost constant within a very short 
time. So we assume C E {1, 1) — Co for all time. (2) C E {1, 1') 
exponentially decays with the time-interval |/ — 1'\ with an 
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exponent 7, i.e., C#(/,/') — Coexp(— -y | Z — l'\). Substitut- 
ing this in Eq. (5), and performing the summations we 
get [13]: 



f(n) * ~C 



n coth 



(?) 



1 — exp(— 771) 



sinh 7 — 1 

Hence, the rate of change of f(n) with time n is 
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7 exp(— 771) 
sinh 7 — 1 



At the beginning, when n is very small, 
A/(n) 1 / C o7 2 



2 \ sinh 7 — 1 



n + const. 



(7) 



(8) 



(9) 



i.e., f(n) evolves quadratically with time (f(n) oc n 2 ). On 
the other hand, when n is very large, 
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This suggests long time linear behavior of f(n). Here we 
identify two distinct time-dependent regimes of the deco- 
herence function f(n) : short time quadratic evolution, 
and long time linear evolution. Identical behavior of the 
decoherence function is reported for the the environment 
consisting of a collection of harmonic oscillators [10]. Thus 
we finally arrive at our goal, and prove that quantum 
chaotic system and a collection of harmonic oscillators be- 
have equivalently as a model of decohering environment. 

We now substantiate our results with numerics. Here we 
consider chaotic kicked top as the model of environment. 
The total Hamiltonian is H(t) = H](t) + Hsif) where 
Hi(t) is the system-environment interaction Hamiltonian, 
and He(£) is the kicked top Hamiltonian. As the environ- 
ment, following version of the kicked top model is used: 



H E (t) = (3 J; 



(11) 



where Ji is the i-th component of the angular momentum 
operator of the top, j is the size of the spin (here, j — 100), 
k is the parameter which decides chaoticity in the system. 
The most popular version of the kicked top model does 
not contain the /3-term. We introduce this extra term 
to remove the parity symmetry i?iJ(i)i? _1 = Hit), where 
R = expiiirJy). The Hilbert space dimension of the kicked 
top is N = 2j + 1 = 201. We set the chaotic parameter at 
k = 99.0 which corresponds to classically strongly chaotic 
system, and the parameter (3 at 0.47. The interaction 
Hamiltonian between the two qubits and the kicked top 
is: 

H T (t) = ah a ® J z . (12) 

The parameter a determines the system-environment 
coupling strength. The system consists of two non- 
interacting qubits, and its coupling agent is h s = S 
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Fig. 1: (Color online) Evolution of the system purity P(n), 
and the parameter A determines the entanglement in the sys- 
tem are plotted for three different system-environment coupling 
strengths a = 0.0001, 0.001, 0.005. Here the initial state is the 
Bell state (see the text). The dotted line is the exact numerical 
calculation, whereas the solid line is the theoretical estimation. 
We estimate P(n) and A from the approximate formula of the 
RDM ps of the system given in Eq. (4). 



j(2) + /(l) <g) sf\ where S ( z l) = /2 and of is third 
Pauli matrix. Hence, the computational basis states 
{|00), 1 01 } , 1 10} , |11)} are also the eigenbasis of h s with 
eigenvalues {—1,0,0,1}. The states 1 01} and 1 10} are the 
degenerate eigenstates of h s . Note that, according to Eq. 
(4), the off-diagonal terms of the system state, correspond- 
ing to this degenerate subspace, will not decay with time. 

Decoherence of the system is studied for two different 
initial states : (1) a Bell state \(f> s ) = 4|(|00) - |11», and 

(2) a product state |^.) = ^(|0)-|1))(8^(|0) + |1)). For 
all numerical studies, a generalized SU(2) coherent state 
[11] is used as initial state for the kicked top environment. 

Fig. 1 presents the result for the Bell state case. In 
the computational basis, diagonal elements of the ini- 
tial density matrix ps(0) are (1/2,0,0, 1/2), and two off- 
diagonal elements have non-zero values : (00|ps(0)|ll) = 
(ll|ps(0)|00) = —1/2. These off-diagonal elements are 
not the part of the degenerate subspace, and therefore 
they decay in time. So in the asymptotic limit : pg ~ 
diag(l/2,0,0,l/2), the purity P and A of this RDM are 
1/2 and 0, respectively. In Fig. 1, we see that the purity P 
and A have reached very close to these asymptotic values. 
The rate, at which these two quantities have reached the 
asymptotic values, is determined by the coupling strength 
a. As expected, for weaker coupling the rate is slower 
than stronger coupling. The most important fact is that 
our approximate formula for the RDM derived from the 
perturbation theory, is not only working well for the weak 
coupling case but also agrees very well for the stronger 
coupling cases. 

Fig. 2 presents the result for the product state case. Ini- 
tially A = 0, and then it becomes negative till n — 1142 
for a = 0.0001, n = 81 for a = 0.001, and n = 4 
for a = 0.005. Upto this time, entanglement between 
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Fig. 2: (Color online) Evolution of P(n) and A are plotted. 
Here initial state is the product state. Other parameters and 
descriptions are same as in Fig. 1. 

two qubits determined by the concurrence was zero. Af- 
ter that, there is a time interval when A > 0, and the 
two qubits become entangled due to their common inter- 
action with the chaotic environment. This shows that, 
like many-body environment [18], quantum chaotic en- 
vironment can also generate entanglement between two 
non-interacting qubits which have no entanglement ini- 
tially. In the asymptotic limit, diagonal elements of ps 
: (1/4,1/4,-1/4,-1.4), and the off-diagonal elements : 
(01|p s |10) = (10|p s |01) = -1/4 corresponding to the de- 
generate subspace, survive from decoherence. Therefore, 
in this limit, the purity P = 3/8 = 0.375, and A = 0. Fig- 
ure shows again a good agreement between the numerics 
and the theory. 

Our theoretical formulation establishes the equivalence 
between quantum chaotic environment and many-body 
environment by showing the equivalence between their 
decoherence function f(n). Therefore, it is very impor- 
tant to investigate f(n) of the kicked top carefully. In 
Fig. 3, we plot the evolution of the decoherence func- 
tion of the kicked top. Our numerical experiment sug- 
gests that the phenomenological formula for f(n) given 
in Eq. (7) is not an exact formula for all time n. It 
describes short time quadratic and long time linear be- 
havior well. So instead of fitting the numerics by the 
phenomenological formula, we interpolate it with a sim- 
ple function f(n) = an + bn 2 , a combination of a linear 
and a quadratic term. The interpolation (solid line of Fig. 
3) gives a ~ 0.1776 and b ~ 1.1561 x 10~ 3 . We expect 
the parameter a to be equal to the coefficient of the lin- 
ear term of Eq. (7), i.e. ^Co coth(7/2). For the chaotic 

kicked top, C = (J z (l) 2 )-(j z (l)) 2 ~ 1/3. In addition, in 
strong chaos regime, the parameter 7 is very large which 
leads to coth(7/2) ~ 1. Hence, the coefficient of the lin- 
ear term of Eq. (7) is approximately equal to 1/6, which 
is very close to the value obtained for the parameter a, 
i.e., a ~ 0.1776. Thus we show that the decoherence func- 
tion of the chaotic kicked top has the properties which are 
similar to the properties of any other model of decohering 
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Fig. 3: (Color online) Evolution of the decoherence function 
f(n) of the kicked top. The dotted line is the exact numerics, 
and the solid line is the theoretical estimation. Inset is showing 
the same in semi-logarithmic scale, where the time-axis is in 
logarithmic scale. 

environment. 

In conclusion, as a model of decohering environment, we 
establish the equivalence between quantum chaotic system 
and any other many-body systems. This is an important 
result, which suggests that instead of complex many-body 
system one can use a simple chaotic system as a model 
of environment. Consequently, analytical investigations 
of environment induced decoherence would be easier, and 
from numerical point of view as well, this simple model 
requires very less computational resources. We substanti- 
ate our analytical formulation with strong numerical evi- 
dences. These show that quantum chaotic system is good 
decoherer, and it can also create entanglement between 
two non-interacting particles. The later result is very im- 
portant due to the recent identification of entanglement as 
a resource for quantum computational and informational 
protocols [2]. 

We thank Drs. A. Tanaka, D. Cohen, and A. Lakshmi- 
narayan for useful comments and discussions. 
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